A Stochastic Approach for Simulating Spatially Inhomogeneous Coagulation Dynamics in the Gelation Regime

We present a stochastic approach for the simulation of coagulation-diffusion dynamics in the gelation regime. The method couples the mass flow algorithm for coagulation processes with a stochastic variant of the diffusion-velocity method in a discretized framework. The simulation of the stochastic processes occurs according to an optimized implementation of the principle of grouping the possible events. A full simulation of a particle system driven by coagulation-diffusion dynamics is performed with a high degree of accuracy, which allows a qualitative and quantitative analysis of the behaviour of the system. Its performance becomes more evident especially in the gelation regime, where the computations become usually very time-consuming.